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Abstract 

We created 4 p-Si models and 4 p-Si:H models all with 50% poros- 
ity. The models contain 32, 108, 256 and 500 silicon atoms with a pore 
parallel to one of the simulational cell axes and a regular cross-section. 
We obtained the densities of states of our models by means of ab ini- 
tio computational methods. We wrote a code to simulate the emission 
spectra of our structures considering particular excitations an decay 
conditions. After comparing the simulated spectra with the experi- 
mental results, we observe that the position of the maximum of the 
emission spectra might be related with the size of the silicon backbone 
for the p-Si models as the quantum confinement models say and with 
the hydrogen concentration for the p-Si:H structures. We conclude 
that the quantum confinement model can be used to explain the emis- 
sion of the p-Si structures but, in the case of the p-Si:H models it is 
necessary to consider others theories. 
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1 Introduction 



Since the discovery of the p-Si photoluminiscence (PL) at room temperature 
reahzed by L.T. Canham in 1990 [Ij, a lot of studies have been carried out 
to understand the origin of this property. The possibihty that p-Si could 
be applied to photonic and optoelectronic devices, motivated the study of 
the dependency of its PL to the temperature, chemical attack variables and 
molecules on its surface. Despite the big amount of scientific reports about 
the p-Si PL, until now, we have not a model that reproduce the characteris- 
tics of all experimental PL spectra reported in the bibliography. 

During the research of the origins of the p-Si PL, some models about the 
nature of the PL were proposed; CuUis et al\2\ say that all these models 
could be reorganized in 6 mayor groups. We only review the fundamentals 
of the confinement quantum model (QC), and of the hydrogen-groups model 
(HS) because we will use them to interpret our results. 

The QC model proposed by Canham in 1990 was the first model created 
with the aim to explain the efficient PL of the p-Si at room temperature. 
In this model, because of the effects of the quantum confinement, the en- 
ergy gap of the silicon nanostructures increase, and their photoluminescence 
appear in the visible range of the electromagnetic spectrum. In accord to 
the QC model, if the size of the silicon nanostructures (where the recombi- 
nation process take place) decrease, then the energy of the photons created 
increase. The QC model is the most successful theory about the nature of 
the p-Si photoluminescence, because it explain various experimental obser- 
vations, with an exception: The PL spectra of oxidized p-Si structures. 

In the other hand, the HS model consider that the hydrogen atoms on 
the p-Si surface (which satisfy the dangling bonds) , could be involved in the 
origin of PL of this material. The HS model was proposed after observing 
that if the hydrogen atoms are removed from the p-Si surface, by a thermic 
procedure, then the intensity of the PL decrease dramatically (but still exist 
because of the QC model argument). Since the p-Si with a rusty surface 
present an intense PL, we understand that, the porous silicon with an hydro- 
gen passivated surface {p-Si:H) is just one of the luminescent forms of the 
p-Si (Cullis et al). 
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Photoluminescence is a difficult phenomenon to simulate, since the dy- 
namics of the excitation and decay of the electrons require precise calculations 
of the excited states and decaying mechanisms. Also, metastable states may 
take part in the process with short, but relevant, mean free lives that require 
time-dependent quantum mechanical calculations [3], [7]. Instead, we have 
decided to initially investigate the photoemission spectrum as a crude ap- 
proximation to the photoluminescence phenomenon. 

Bearing in mind the QC and HS models, we propose a simple photoe- 
mission model in which the excited electrons of the conduction band, decay 
successively into the unoccupied electronic states of the valence band (see 
figure HI). The emission spectra that we obtain for the p-Si and p-Si:H struc- 
tures are comparable to the experimental spectra reported by Cullis et al. 

We are aware that our simulated p-Si structures do not have experimental 
counterpart, because the number of silicon atoms involved in each simulation 
is finite due to the computational cost of the ab initio calculations. Also, in 
our simulations is not possible to reproduce all the chemical molecules that 
generate after the chemical attack of the p-Si surface, because depending of 
the atmosphere where the attack was realized, molecules or atoms besides 
hydrogen atoms, could satisfy the dangling bonds of the p-Si surface modi- 
fying the form of its PL spectrum. 

However, until now the doubts about the origins of the p-Si PL remains, 
so, our work contributes to calculate the photoemission spectra of p-Si and 
p-Si:H structures with ab initio molecular dynamics methods. In the follow- 
ing section, we describe the construction method used to generate the p-Si 
and p-Si:H structures, and the way we obtain their photoemission spectra. 



2 Method 

The p-Si structures were created using the Cerius2 interface [8j. We repro- 
duced the silicon cell diamond-like which contains 8 silicon atoms 2 times 
along each of the coordinate axes to obtain a 64 silicon-atom supercell, 3 
times to increase the number of atoms to 216, 4 times for a supercell with 
512 atoms and finally, 5 times along each direction to obtain 1000 atoms in 
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the silicon supercell. 



Once the superceUs were created, we ehminated 50% of their sihcon atoms 
carving a central pore parallel to the z-axis and with a regular cross section. 
After eliminating 50% of the silicon atoms contained in each supercell, we 
obtained 4 p-Si structures with 32, 108, 256 and 500 silicon atoms. The pore 
of the p-Si structure with 32 silicon atoms has an irregular transversal section 
because the great number of atoms removed from the supercell. In the other 
hand, the pore of the supercell with 500 silicon atoms shows a regular cross 
section because even if we remove 50% of their atoms, its silicon backbone is 
bigger than the one for the 32 silicon atoms model. 

The 4 structures of p-Si were subjected to a geometry optimization pro- 
cess with the objective of finding their equilibrium configurations (see figures 
l(a)| to l(d)[ ). We perform the geometry optimization process with FAST 



STRUCTURE a code included in the CERIUS2 interface based in the Lin 
and Harris approach [5j , j4j , and the VWN [6j functional of correlation energy. 

In table [H we indicate the computational parameters used for the geom- 
etry optimizations of the 4 p-Si structures. The integration grid refers to 
the space were the atomic functions, atomic densities. Coulomb potential, 
exchange and correlation potentials are represented. To calculate the hamil- 
tonian matrix elements it is necessary to put each atom in the center of the 
integration grid. A coarse integration grid leads to a poor representation 
of the atomic functions, whereas a fine grid is the best approximation but 
its computational cost is elevated. The minimal bases indicate that we use 
the valence atomic orbital to construct the wave function, and the minimal 
density indicates that the electronic density should be obtained just using 
the valence atomic orbital as well. Finally, the cut-off parameter corresponds 
to the distance at which the wave function vanish. 

At the end of the geometry optimizations of our p-Si structures, the in- 
terface brings forth their radial distribution functions (figures 2(a) to |2(d) 



and [3]) and densities of states. We used part of the values of the densities 
of states of the p-Si models as the input information for our emission code 
(EDOxEDO.f). We consider that the electronic states of the conduction den- 
sity are unoccupied and those of the valence density are occupied. Following 
the decay order proposed by our code (figured]), the excited electrons of the 
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Parameter 


Option 


Integration grid 


Coarse 


Bases 


Minimal 


Density 


Minimal 


Cut off radius 


5 A 



Table 1: . Computational parameters used for the geometry optimization 
process implemented by the FAST STRUCTURE code for the p-Si and 
p-Si:H models. 



conduction states decay into the valence density states emitting a photon, 
whose energy is equal to the energetic difference between the two electronic 
states involved. The energies of the photons are computed as an histogram 
and then, we adjust a gaussian curve to the histogram to compare it with the 
p-Si PL experimental spectrum reported by CuUis et al (CuUis et al). To 
calculate the emission of all our structures, we only consider the electronic 
states that belong to a 2 eV interval in the valence and conduction densities 
of states. We decide to work out with a 2 eV interval states because ex- 
perimentally, the majority of the PL excitation sources have energies of this 
order of magnitude. 



The p-Si:H models were constructed taking the p-Si optimized structures 
as the initial atomic arrangement. We saturated with hydrogen the p-Si opti- 
mized structures. After the saturation with hydrogen, we applied a geometry 
optimization to the p-Si:H structures using the parameters of table [H In fig- 
ures 5(a) to 5(d) I we present the 4 optimized structures of p-Si:H. We obtain 
the valence and conduction states from the densities of states of the p-Si:H 
models in the same way than for the p-Si structures. Again, to obtain the 
emission spectra of the p-Si:H models we just consider the electronic states, 
in the valence and conduction densities, that are in an 2 eV interval. The 
order of decay of the excited electronic states into the unoccupied valence 
states is shown in figure [H The distribution of the photon energies is ob- 
tained after a gaussian fit. 



In the results section we exhibit the photoemission spectra of the p-Si 
and p-Si:H models. We compare our photoemission spectra with the PL 
spectrum reported by CuUis et al. 
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(a) FDRpSi32 




(b) FDRpSilOS 




1D r= 1-4 ia 




(c) FDRpSi256 



(d) FDRpSiSOO 



Figure 2: Radial distribution functions of the p-Si optimized models 
with 32 (a), 108 (b), 256 (c) and 500 (d) silicon atoms. 
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Figure 3: . Interatomic distance variation of the p-Si optimized models, 
as a function of the size of their sihcon backbones. 
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Figure 4: . Order (considered in the EDOxEDO.f code.), of the excited 
electrons decay from the conduction band to the valence band. 
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(a) pSi32H42 (b) pSil08H68 




(c) pSi256H153 (d) pSi500H196 

Figure 5: . p-Si:H optimized structures with 32 silicon atoms and 42 of 
hydrogen (a), 108 sihcon atoms and 68 of hydrogen (b), 256 sihcon atoms 
and 153 of hydrogen (c) and 500 sihcon atoms and 196 of hydrogen (d). 
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3 Results and discussion 



As we can see from figures 2(a)| to 2(d)[ the radial distribution functions of 



the p-Si structures become similar to the radial distribution function of a 
crystalline silicon structure, until the number of silicon atoms reach 500 (fig- 
ure 2(d)D . The variation of the interatomic distance of the p-Si structures as 
a function of their silicon backbones size is shown in figure [3l The backbone 
of the p-Si structure with 500 silicon atoms is thicker than the one of the 
model with just 32 silicon atoms, therefore, the thinner silicon backbone is 
really deformed after a geometry optimization process because of the pres- 
ence of the porous in the simulation cell, then the most probable distance 
between silicon neighbors increase to 2.55 A. In the other hand, the backbone 
of the p-Si structure with 500 silicon atoms, is thick enough to remain as a 
crystalline silicon structure despite the central porous, so then the distance 
between first neighbors is 2.35 A, as for a crystalline silicon structure. 



Figures [6(a) , 6(b), 6(c) and |6(d) show the gaussian curves fitted to the 



histogram of each p-Si model obtained by our EDOxEDO.f code. In figure 
7(a)| we present the emission spectra (just the gaussian curves) of the 4 p-Si 



optimized structures; for these 4 structures their photoemission energy range 
from 1.2 eV to 2.8 eV, its mean maximum photoemission energy is 1.94 eV, 
and as we can observe, the intensity of the emission increase if the size of the 
backbone increase, but its maximum energy decrease in accord to the results 
of Cullis et al. 



In figure 7(b)| we compare the p-Si PL spectrum reported by Cullis et al 



and the emission spectrum of our p-Si model with 500 silicon atoms; from this 
figure we can appreciate that our spectrum is larger than the Cullis et al but 
both have a similar form. Also, the maximum photoemission energy of our 
spectrum is 0.34 eV to the blue of the electromagnetic spectrum considering 
the va 

In figure [8] we superpose the emission spectra of the 4 models of p-Si:H 
calculated by the EDOxEDO.f code. We did not find experimental photoe- 
mission spectra for the p-Si:H to compare to our results. 
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(c) pSi256 (d) pSiSOO 

Figure 6: Gaussian fitting of the p-Si model histograms (a) with 32, (b) 
108, (c) 256, and (d) 500 sihcon atoms. 




(a) Emission spectra of the p-Si models (b) Simulation and experiment 

Figure 7: (a) Emission spectra of the 4 p-Si structures calculated with 
the EDOxEDO.f code, (b) p-Si PL spectrum reported by Cullis et al. 
(black curve) and the emission spectrum of our p-Si model with 500 silicon 
atoms (red curve). 



11 




LiiEiH-.! ! ty "I 

Figure 8: . Emission spectra of the four p-Si:H models calculated with 
de EDOxEDO.f code. 



4 Conclusions 

We observed that the form of the emission spectra, obtained after a gaussian 
fitting of the histograms calculated by the EDOxEDO.f code, is similar to the 
form of the CuUis et al spectrum. The maximum width of our spectra is 1.6 
eV, their minimal energy is 1.2 eV and the maximal one is 2.8 eV. The mean 
emission energy calculated was 1.94 eV, this is, 0.34 eV toward the blue of 
the electromagnetic spectrum, with respect the CuUis et ai maximum energy. 



Is clear from the figure [7(a) that decreasing the number of silicon atoms in 



the simulation cell, moves toward bigger energies the maximal energy of the 
emission spectra. We think that this fact is due to the size of the quantum 
confinement region in the silicon backbone; for example, in the 32 p-Si struc- 
ture, the silicon backbone is smaller than the one of the 500 silicon atoms 
model and, in accord with the quantum confinement model, his emission en- 
ergy should be bigger than the one of the p-Si 500 atoms model. 



From figure 7(a) we conclude that decreasing the number of silicon atoms 
in the p-Si supercells results in a smaller emission intensity. We think that 
because the p-Si 500 silicon atoms model has a crystalline structure ( |2(d) ) 



compared with the 32 atoms structure ( 2(a)D , the radiative recombination 



of the charges dominates in the crystalline structures of the p-Si and is not 
significant in the amorphous ones. 
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We observe from figure 7(b)| that the emission spectra of our p-Si 500 sil- 
icon atoms model is toward the right of the electromagnetic spectrum, with 
respect to the PL maximum reported by Cullis et al. We consider that the 
size of the supercell is the direct cause and probably, if we could be able 
to increase the number of atoms, to thousand of they to construct the p-Si 
models, theirs emission maximus will be near the experimental ones reported. 
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